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Abstract 

Wc describe how density-functional theory, well-known for its many 
uses in ab initio calculations of electronic structure, can be used to 
study the ground state of inhomogeneous model Hamiltonians. The 
basic ideas and concepts are discussed for the particular case of the 
Heisenberg model. As representative applications, illustrating scope 
and limitations of the procedure, we calculate the ground-state energy 
of one-, two- and three-dimensional antifcrromagnetic Heisenberg mod- 
els in the presence of boundaries and of impurities in the bulk and at 
the surfaces. Correlations are shown to lift degeneracies present in the 
mean-field approximation. Comparison with exact (brute force) diag- 
onalization shows that the density-functional results are a significant 
improvement over the mean-field ones, at negligible extra computa- 
tional cost. 



"This is the (slightly modified) English translation of a chapter, originally written in 
Portuguese, of a proceedings volume on density-functional theory: 40 anos DFT: Uma 
perspectiva Brasileira, eds. K. Capelle, A. J. R. da Silva and A. Fazzio. 
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1 Introduction 



Traditional density-functional theory (DFT) is based on the Hohenberg- 
Kohn theorem, according to which all observables of a quantum many- 
body system are determined uniquely by the ground-state charge density 
[1, 2]. These observables are thus functionals of the density, and a ma- 
jor task of DFT is the construction of ever better functionals for calculat- 
ing, e.g., ground-state energies in terms of the ground-state density. The 
local-density approximation (LDA) is one such functional, expressing the 
exchange-correlation (xc) energy of an inhomogeneous many-particle sys- 
tem in terms of its charge density. The density itself is normally calculated 
by means of a self-consistent field (SCF) procedure, solving effective single- 
particle equations known as the Kohn-Sham equations [3]. 

While widely used in electronic-structure calculations, both in solid-state 
physics and in quantum chemistry, the concepts and tools of DFT, such as 
the Hohenberg-Kohn theorem, the Kohn-Sham equations, and the LDA and 
its many improvements, are much more general, and can be applied to many 
quantum systems not described by the ab initio (Coulomb) Hamiltonian 
underlying applications of DFT in band-structure and quantum-chemical 
calculations. 

In particular, they can also be applied to model Hamiltonians. In a 
model Hamiltonian description of a physical phenomenon one does not pri- 
marily focus on material-specific details, but on generic features of classes of 
systems. The use of such models has a long tradition in many-body physics 
and statistical mechanics, and the list of common models is long and includes 
famous examples such as the Ising, Heisenberg and Hubbard models. 

DFT can be a useful tool in the study of such models in the presence 
of inhomogeneities [4, 5, 6]. Below we briefly describe the construction and 
some representative applications of DFT for the antiferromagnetic Heisen- 
berg model [5, 7, 8]. To prepare the ground, we present in Sec. 2 a system- 
atic comparison between traditional ab initio (Coulomb) DFT and DFT for 
model Hamiltonians. Section 3 introduces the Heisenberg model. In Sec. 4 
we describe expressions for the ground-state and correlation energy of ho- 
mogeneous (spatially uniform) Heisenberg models. In Sec. 5 we use these 
expressions, in the spirit of an LDA, to construct a local-spin approximation 
for inhomogeneous one, two and three-dimensional Heisenberg models. In 
Sec. 6 these functionals are used to calculate the ground-state energy of in- 
homogeneous chains, ladders and cubes, and compared with mean-field and 
exact diagonalization data. Sec. 7 contains a brief summary and assessment 
of the approach. 
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2 Comparison between ab initio DFT and DFT 
for model Hamiltonians 



In Table 1 we present a schematic comparison between traditional, ab ini- 
tio, DFT and DFT for two widely studied model Hamiltonians denned on 
lattices, the Hubbard and the Heisenberg model. The basic idea of lattice 
DFT for model Hamiltonians goes back to Gunnarsson and Schonhammer 
[9], in the context of the Hubbard model. A viable parametrization for the 
model's xc energy was proposed in [4], while corresponding work for the 
Heisenberg model is based on [5]. The following paragraphs describe Table 
1, line by line. 

Starting point is, in all Hohenberg-Kohn (HK) theorem, estab- 

lishing that an intensive density-like quantity, coupled to a suitable gener- 
ating field, determines the ground-state wave function of the model under 
study. While in ab initio DFT this intensive variable is the charge density 
n(r), in lattice DFT for the Hubbard model it is the site occupation number 
rii, and for the Heisenberg model it is the local spin Sj. 

Having proven a Hohenberg-Kohn theorem, the variational principle al- 
lows us to obtain the ground-state energy and the ground-state density 
distribution by minimizing a density functional. In ab initio (Coulomb) 
DFT and Hubbard DFT this minimization is normally done indirectly, via 
self-consistent solution of Kohn-Sham (KS) equations. This indirect mini- 
mization allows one to circumvent having to approximate the single-particle 
kinetic energy as a density functional. The Heisenberg model does not have 
a kinetic-energy term, and the minimization is most conveniently performed 
directly. 

The xc energy is also unknown. Starting point for most approximation 
schemes for the xc functional is the local-density approximation, in which 
one first defines an auxiliary spatially uniform many-body problem, which 
can be solved more easily. In ab initio DFT this auxiliary many-body system 
is the homogeneous electron liquid, while for lattice DFTs spatial uniformity 
means equivalent sites. The solution of this uniform auxiliary problem may 
be very difficult, but it is still a lot simpler than the solution of the inho- 
mogeneous many-body problem. In ab initio DFT the best values for the 
xc energy come from Quantum Monte Carlo simulations [10]. The one- 
dimensional Hubbard model, and the one-dimensional spin one-half Heisen- 
berg model, have formally exact solutions by means of the Bethe Ansatz. 
For these systems the xc energy can be extracted, in principle, exactly. In 
higher dimensions or for larger spins this can only be done approximately. 

Having solved the spatially uniform system, the usual next step is to 
parametrize the xc energy in a form comvenient for numerical calculations. 
In ab initio calculations common parametrizations are those of Vosko, Wilk 
and Nusair (VWN) [11], Perdew and Zunger (PZ) [12], and Perdew and 
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Table 1: Diferent realizations of the concepts of DFT. See the main text 
for a description of the various entries. 
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Wang (PW) [13]. For the Hubbard model a viable parametrization was 
proposed in [4]. Alternative approaches are discussed in [9, 14]. For the 
Heisenberg model several parametrizations were proposed in [5]. Some of 
these are briefly described below. 

Once a parametrization of the per-site xc energy of the uniform model is 
available, it can be applied site by site in the corresponding inhomogeneous 
model. In ab initio calculations the inhomogeneity arises from the crystal 
lattice in a solid, or the atomic nuclei in a molecule. In lattice models it 
arises from inequivalent sites. For the Heisenberg model three examples of 
such lattice-DFT calculations for models with inequivalent sites are given 
below. These examples complement those presented in [7]. 



3 The Heisenberg model 

One of the requirements on a formalism purporting to describe the magnetic 
properties of matter is that it should predict the existence of long-range 
order of ferro or antiferromagnetic nature [15, 16]. Within DFT, this re- 
quirement is met by spin-density- functional theory (SDFT). Within model 
Hamiltonians, the simplest model displaying such order is the Heisenberg 
model. 

(ij) 
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in which is the " exchange interaction" between spins at nearest-neighbour 
sites. The terminology "exchange" here arises from the original motiva- 
tion for writing Eq. (1) as an effective model accounting for the conse- 
quences of the Pauli exclusion principle [17, 18]. However, if J« is cal- 
culated from many-body theory or fitted to experiment, it generally also 
contains contributions from correlation, i.e., many-body effects beyond the 
single-determinant (Hartree-Fock like) approximation. 

Depending on the sign of Jij and its dependence on the site labels i and j, 
the model (1) can have ferromagnetic, antiferromagnetic or noncollinear spin 
configurations [16]. Below we focus on the antiferromagnetic case, which is 
relevant for a large class of magnetic insulators, among which the undoped 
state of cuprate superconductors may be one of the most important [19]. 
However, the applications of model (1) transcend the boundaries of physics, 
as illustrated by chemical applications to the magnetic states of conjugated 
hydrocarbons and organometallic compounds [19-24]. 

Although almost 80 years have passed since the model (1) was first pro- 
posed, a complete solution is known only for the limiting cases of few sites 
and, by means of the Bethe Ansatz (BA), for spin S = 1/2 in one dimen- 
sion [26, 27]. Complete numerical diagonalization is possible for up to a 
few dozens of sites. The density-matrix renormalization group (DMRG) 
and Quantum Monte Carlo (QMC) techniques can be applied also to larger 
systems, but at considerable computational cost. 

Moreover, all these techniques, BA, DMRG, QMC, etc. encounter signif- 
icant difficulties in spatially nonuniform systems, with broken translational 
symmetry, such as in the presence of impurities, boundaries or externally 
applied nonuniform fields. A modification of the homogeneous Heisenberg 
model describing this type of inhomogeneity is 



where Bj is an external or internal magnetic field, which can vary from one 
site to the next, and Jy extends the spin-spin coupling to other than nearest 
neighbours. 

In these situations DFT provides a computationally viable approach to 
the ground-state energy and density, and their dependence on the system 
parameters. However, as explained in the preceding section, already the 
simplest practical approximation to DFT, the LDA, requires as an input 
the solution of the uniform system. For this reason we turn, in the next 
section, to a brief discussion of the ground state and correlation energy of 
the homogeneous Heisenberg model. 




(2) 
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4 Ground-state energy of the homogeneous Heisen- 
berg model 

The simplest approximation to Hamiltonian (1) is the mean-field (MF) one, 
in which the vector operators S are substituted by classical (commuting) 
vectors S. In this approximation the per-spin ground-state energy of a 
homogeneous linear, quadratic or cubic lattice comprising N sites with spin 
S, in d dimensions, is 

e M F{s) = Eig& = _ JdS 2 (3) 

An improved estimate of this energy can be obtained from spin-wave (SW) 
theory [28], according to which 

ejWfS) = -Jd S 2 + Jd'^H- - 1)5 = e% IF (S) + Jd~ l l\- - 1)5 . (4) 

7T 7T 

(Here we have used the conjecture of Ref. [5] to write the last term as a closed 
function of d.) As an example, for d = 1 and 5 = 1/2, spin-wave theory 
predicts e^ w = —0.431690 J, which is about 2.6% off the exact Bethe-Ansatz 
result efi A = (1/4 -In 2) J. 

The correlation energy e c (S) is defined as the difference eo(S) — e^f F (S), 
where eo(5) is the exact ground-state energy. Within SW theory, we have 
thus 

e s c W(S) = e s w (S) - e™ F (S) = Jd~ l l\- - 1)S . (5) 

7T 

A more precise expression for e c (S) in d = 1 has recently been proposed 
[5, 29] on the basis of previously obtained density- matrix renormalization 
group data [29], 1 

e° MRG (S) e s c w {S) 0.0030 0.0015 

^-7^ = ~ °- 03262 " — + — 

0.338 - ^ + ^) e^ s cos(2vr5) . (6) 

D b J 

For 5 = 1/2 this expression predicts e$ = —0.446253 J, which deviates 
only by 0.7% from the exact result. 

The utility of the MF, SW and DMRG expressions for eo(S) is, a priori, 
rather limited: MF theory is not reliable due to its neglect of correlation, 
and the SW and DMRG expressions are only applicable to spatially uniform 
systems, in which b is a constant. (The DMRG expression is also limited 
to d = 1.) Restriction to spatial homogeneity is a serious limitation, since 

The fit originally proposed in Ref. [29] did not include the cubic terms in 1/5 and 
did not recover the Bethe-Ansatz exact result at S = 1/2. The cubic terms were added in 
Ref. [5] in order to extend the fit to S = 1/2. 
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spatial inhomogeneities are ubiquitous in real systems. Examples are effects 
of external magnetic fields, magneto-crystalline anisotropy, boundaries, im- 
purities, defects etc. Such situations are hard to deal with by traditional 
methods because translational symmetry is broken. To overcome this diffi- 
culty we employ density- functional theory [2, 30], within which expressions 
for the uniform system find applications as input for LDA-type approxima- 
tions for inhomogeneous systems. 

5 Ground-state energy of the inhomogeneous Heisen- 
berg model: spin-distribution functionals and 
local-spin approximation 

The Hohenberg-Kohn theorem for the Heisenberg model, proved in Ref. [5], 
shows that the expectation value of any observable O of that model is a 
functional of the ground-state spin distribution Sj, given, in principle, as 
OfSj] = (^[SjJIOI^fSj]), where Sj = ( x I / |Sj| x I / ) is a classical spin vector and 
not an operator. This functional is universal with respect to external fields 
Bj, but not with respect to changes in the spin-spin interaction Jij. 

For the Heisenberg model the ground-state energy and its minimizing 
spin distribution are obtained most conveniently by direct minimization of 
the functional 



which for homogeneous antiferromagnetic systems on linear, square and cu- 
bic lattices reduces to Eq. (3). To approximate £7 c [Sj] we employ the LDA 
concept, applied directly to the spin vectors. The resulting local-spin ap- 
proximation (LSA) consists in substituting, locally, the variable S in the 
expressions for the homogeneous system, by the variable |Sj|, according to 




(7) 




(8) 



{ij) 



E C (S) 



(9) 



N 



Within the SW approximation for E c we thus have 




(10) 



so that the ground-state energy functional in the hSA sw becomes 



E LSA-sw [Si] =J J2S i -S j + J d~^(- - 1) E is*l • (11) 



{ij) 
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Figure 1: Per-site ground-state energy of an open 10-site chain of spin 1/2 
atoms, obtained exactly, in the mean-field approximation, and with LSA. 
The inset illustrates the difference between LSA SW and LSA DMRG . 

Similarly, the LSA DMRG is 

E LSA- D MRG [Si] = j£ Si .S i + £ e r* G (S)ls-^| , (12) 

(y) (v) 
where e® MRG {S) is defined in Eq. (6). 

6 Applications: Finite-size quantum spin chains, 
ladders and cubes 

As an application of these approximations we now consider 10 spin S sites 
on a chain with open boundary conditions. As a consequence of the bound- 
aries this is already an inhomogeneous system, in which not all sites are 
equivalent. In Fig. 1 we compare values for the per-site ground-state energy 
obtained from the mean-field approximation, with values obtained from the 
LSA'W an( j LgA. £>Mi?G functionals, for some values of the spin S. We also 
include exact values obtained from numerical diagonalization of Hamiltonian 
(1). The comparison of the three approximation schemes with the exact 
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Figure 2: Per-site ground-state energy of an N y x N x lattice of S = 1/2 spins, 
with open boundary conditions. N p = N y x N x is the total number of sites, 
ranging from 2 to 10 4 . Curves and data pertaining to a one-dimensional 
chain of length N x are in blue, those pertaining to two-legged ladders of size 
2 x N x in red, and the curve for multi-legged ladders of size 100 x N x is in 
green. See main text for more detailed explanation of the different data sets 
and curves. 



data illustrates the significant improvement obtained by means of the LSA, 
as compared to the MF approximation. It also illustrates that the LSA ,s,/l/ 
and LSA are very similar on this scale. In the inset we plot the differ- 

ence — e$ w versus 5, to illustrate the scale of the difference between 

both formulations of the LSA. For larger systems than L = 10, and/or larger 
spins than 5 = 2, full numerical diagonalization becomes computationally 
very expensive. On the other hand, cheap mean-field calculations are not 
reliable. The LSA concept is thus a useful tool for estimating energies of 
systems that are too large for exact calculations, at much better accuracy 
than obtained from mean-field calculations. 

The main utility of the LSA concept, however, lies in the ease with which 
it is extended to higher dimensions or to more inhomogeneous systems. As a 
second illustrative application of the LSA we thus present, in Fig. 2, results 
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Table 2: Ground-state energy —Eq/NJ in the mean-field approximation, 
of a 4 x 4 x 4 cube with 63 spin 1/2 sites and one spin Si site located on a 
corner, an edge, a face or in the bulk. 





1/2 


1 


3/2 


2 


5/2 


3 


corner 


0.5625 


0.5742 


0.5859 


0.5977 


0.6094 


0.6211 


edge 


0.5625 


0.5781 


0.5938 


0.6094 


0.6250 


0.6406 


face 


0.5625 


0.5820 


0.6016 


0.6211 


0.6494 


0.6602 


bulk 


0.5625 


0.5859 


0.6094 


0.6328 


0.6562 


0.6797 



Table 3: Ground-state energy -E /NJ in the LSA SW of the same system 
described in Table 2 and Fig. 3. Inclusion of correlation lifts the degenera- 
cies predicted by the mean-field approximation. 





1/2 


1 


3/2 


2 


5/2 


3 


corner 


0.7080 


0.7220 


0.7360 


0.7500 


0.7640 


0.7780 


edge 


0.7080 


0.7259 


0.7438 


0.7617 


0.7796 


0.7975 


face 


0.7080 


0.7298 


0.7516 


0.7734 


0.7952 


0.8170 


bulk 


0.7080 


0.7337 


0.7594 


0.7851 


0.8108 


0.8366 



for systems representative of the crossover from one to two dimensions, such 
as spin ladders. Figure 2 contains two sets of exact data. The first, labelled 
exact rf=1 , is the per-site energy of a linear open chain with N x sites. For 
this same system we also show results from the mean-field approximation, 
Eq. (8) (labelled MF d=1 ), and the LSA DMRG and LSA SW approaches (la- 
belled LSA-DMRG d=1 and LSA-SW d=1 , respectively). The arrow labelled 
0.443 indicates the exact value in the thermodynamic limit N x — > oo, ob- 
tained from the Bethe-Ansatz. Comparing these data we see that the MF 
approximation fails badly, while LSA DMRG provides an excellent approxi- 
mation to the exact data for not too small N. 

The second set of exact data in Fig. 2, labelled exact 15 , represents 
the energy of spin ladders of size 2 x N x , as a function of N x . Again, we 
compare the exact data with results from the mean field (MF LAD ) and both 
local-spin (LSA-DMRG LAD and LSA-SW d=1 ) approximations. The arrow 
at 0.578 is the value extrapolated for N x — > oo from Quantum Monte Carlo 
data [31]. Although for the ladders LSA DMRG and LSA W are not as good 
as for chains (which is easy to understand, considering the origin of their 
correlation energy in a one-dimensional system), it is still a significant im- 
provement on the mean-field approximation, which for the two-legged ladder 
predicts energy values that fall closer to the exact ones for the chain than 
to those for the ladder itself. In numbers, at N x = 12 the MF approxi- 
mation differs by 37% from the exact value, while the ~LSA DMRG value is 
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Figure 3: 4x4x4 cube of spin 1/2 with a substitutional impurity at 
the corner (illustrated), edge, face or in the bulk. The mean-field ground- 
state energies resulting from the different impurity positions show certain 
degeneracies that are lifted by inclusion of correlations. (See data in Table 
2 and 3 and discussion in main text.) 



off by 2.4%. Still in Fig. 2, the curve labelled LSA-SW d=2 represents val- 
ues for more and more two-dimensional systems, of size N x x 100, obtained 
by combining the two-dimensional mean-field approximation with the two- 
dimensional SW approximation as an input for LSA. The value obtained 
from LSA 5M/ at N x = 100 deviates by only 3.2% from the asymptotic result 
0.6696 J, estimated in Ref. [32]. 

As our third example we consider a three-dimensional cubic arrangement 
of spins with 4 x 4 x 4 = 64 spin 1/2 sites, one of which is replaced by a sub- 
stitutional impurity of spin Si. This impurity can sit in the bulk of the cube 
(8 equivalent positions), on one of its faces (24 equivalent positions), on an 
edge (24 equivalent positions), or on the corner sites (8 equivalent positions). 
Fig. 3 illustrates this system, and Tables 2 and 3 contain the corresponding 
mean- field and LSA SW ground-state energies for several values of Si. We 
note that the mean-field approximation predicts certain degeneracies that 
are lifted on including correlation via the LSA: in the MF approximation an 
Si = 2 impurity on a face of the cube and an Si = 3 on a corner site are 
degenerate (both having energy e^f F = — 0.6211J), as are Si = 1 impurity 
in the bulk and an Si = 3/2 one at a corner {e^ F = — 0.5859J). A triple 
degeneracy occurs between an Si = 5/2 impurity at a corner, an Si = 2 one 
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on an edge site, and a Si = 3/2 one in the bulk. (eQ IF = — 0.6094J). All 
of these degeneracies are lifted by correlations, which shows that they are 
artifacts of the mean-field approximation. This type of information could be 
useful for the design of self-assembled magnetic nanostructures, in which the 
magnetic atoms will tend to spontaneously occupy the energetically favoured 
sites. Other applications of the LSA SW and LSA DMRG functionals can be 
found in Refs. [5, 7, 8]. 

7 Summary 

The above examples show that DFT within the local-spin approximation 
can be a useful tool for calculating ground-state energies of spatially inho- 
mogeneous Heisenberg models. The computational effort is much less than 
for a QMC or DMRG calculation, allowing applications to systems in two 
or three dimensions, and with reduced translational symmetry. 

Exact or nearly exact benchmark values can only be obtained for small 
numbers of sites and simple geometries, or in the thermodynamic limit. LSA 
calculations, on the other hand, are not more expensive, computationally, 
than the mean-field ones, and much more accurate, but this advantage is 
partially offset by the fact that, unlike an exact, BA, or QMC calculation, 
LSA does not yield the many-body wave function. LSA is also less accu- 
rate than DMRG, but more easily applicable to systems in two and three 
dimensions. 

Methodologically, DFT for the Heisenberg model differs in three aspects 
from Coulomb (ab initio) DFT [2] and from DFT for the Hubbard model [4]: 
(i) LSA proceeds exclusively in terms of the spin density and makes no use 
of the charge density (differently from the local-spin-density approximation 
of Coulomb DFT, which uses charge and spin densities), (ii) Numerical 
results are obtained by directly minimizing an energy functional, and not 
via self-consistent solution of effective single-particle equations. In this sense 
Heisenberg-model DFT in the LSA represents a realization of orbital-free- 
DFT [33]. (iii) For the Heisenberg model the dependence of the correlation 
functional on dimensionality is known to a good approximation [5]. This 
dependence is unknown in other formulations of DFT. 

This last observation also represents a different way in which DFT for 
model Hamiltonians can be useful: While the above mostly focused on how 
DFT can be usefully employed to obtain information on model Hamiltoni- 
ans, it is also possible to regard these models as theoretical laboratories in 
which concepts and tools of DFT can be investigated in a simpler and better 
controlled environment than in an ab initio calculation. 
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